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Abstract 


The design of heat exchangers traditionally focuses on the known con- 
straints of the problem such as inlet and outlet temperatures, flow rates, and 
pressure drops. This leads mainly to a sizing problem where the designer 
must select surfaces, flow configuration, and materials to meet the mini- 
mum design objectives. An alternate approach based on an acceptable level 
of thermodynamic irreversibility (entropy generation) has been proposed. 
When the entropy generation level has been set, the geometric parameters 
of the heat exchanger can be determined. The design of a plate-fin type, 
gas-to-gas recuperator for a regenerative open Brayton cycle has been used 
as a demonstrative device. The resulting heat exchanger designs are then 
examined to determine what caused the differences and why either method 


should be preferred over the other. 
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Chapter 1 


INTRODUCTION 


1.1 Background 


Energy conservation is a topic that has recently received considerable atten- 
tion. Available work has been recognized as a valuable commodity and its 
destruction in seemingly otherwise efficient engineering processes has also 
seen a renewed increase in interest. The use of second law analysis and ther- 
modynamic irreversibility minimization have been proposed as techniques 
that should become an integral part of the design of engineering processes 
and components [1]. It is not apparent what the relationship of traditional 


methods of design and these irreversibility minimization techniques will be. 


1.2 Conventional Design Procedures 


The design of a component, such as a heat exchanger, generally involves the 
sizing of that component to meet specified performance parameters within 
known constraints. The heat exchanger area and volume are usually de- 
signed to be the minimum required to meet the specifications as this will 
also usually be the best design from an ecomonic aspect. The fluid inlet and 
outlet temperatures, flow rates, and pressure drops are usually specified and 
it is the designer’s task to determine construction type, flow arrangement, 
materials and surfaces to meet those requirements [2]. The effect of the 


component design on the overall system performance should also be consid- 
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ered to ensure the tradeoffs performed within its design are still valid when 


integrated into the system. 


1.3 Irreversibility Minimization 


Any heat transfer process is generally accompanied by thermodynamic 1r- 
reversibility or entropy generation. The entropy generated is in direct pro- 
portion to the amount of useful work dissipated in that process. If this 
irreversibility can be minimized within a particular component of a power 
cycle, the useful power output of the cycle should increase [3]. Bejan has 
shown how the entropy generation rate can be reduced in a counterflow gas- 
to-gas heat exchanger [4], and it is this method that will be used to create 
designs that will be compared to conventionally designed heat exchangers, 


and their impact on overall cycle performance. 
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Chapter 2 


SELECTION OF 
BRAYTON CYCLE 
PARAMETERS 


2.1 Cycle Components 


A regenerative open Brayton cycle was chosen as the vehicle to test the 
different design methods for the recuperator. A schematic of the cycle com- 


ponents is shown in Figure 2.1. 


Recuperator Combustor Turbine 
l 
Compressor 


Regenerative Cycle Schematic 
Figure 2.1 
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A compressor polytropic efficiency was assumed and used throughout 
all calculations, 7¢-p = 0.95. Similarly, a turbine polytropic efficiency was 
assumed, 7,¢ = 0.90. For ease of calculations, only one fluid was used, 
l.e., air was modeled as a perfect gas with cp = 1.0 kyj/kg —° K and R = 
0.287 kj/kg—° K both assumed constant. The combustor, then, was treated 


as a perfect heat transfer device, but was not an injection point for fuel. 


2.2 Cycle Parameters 


A reasonable compressor inlet temperature was selected, T; = 300°K, along 
with an inlet pressure of P, = 1.126 x 10°N/m?. A nominal value of turbine 
inlet temperature was also selected, Ty, = 1300°K. The selection of an 
appropriate pressure ratio required the following consideration. 

The addition of a recuperator to a simple open Brayton cycle causes a 
shift in the pressure ratio for maximum cycle efficiency from that of the non- 
regenerative cycle. A heat exchanger effectiveness of ¢€ = 0.8 was selected 
and cycle efficiency without losses was calculated for pressure ratios varying 
from 4 to 10. A sample calculation is contained in Appendix A, results are 


plotted in Figure 2.2. 


Cycle Efficiency vs Pressure Ratio 
Figura 2.2 
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Based on these results, a pressure ratio of P2;/P, = 5 was selected. With 
this pressure ratio determined, the maximum cycle efficiency without losses, 
Noo, could be calculated. This quantity will be used for comparison purposes 
and was computed as no, = 0.5760. 


2.3 Heat Exchanger Characteristics 


A plate-fin heat exchanger surface was selected from those cataloged in ref- 
erence [5]. Its characteristics and dimensions are shown in Figure 2.3. It has 
an uncomplicated geometry and a reasonably wide range of Reynolds num- 
bers in which consistent results can be expected. The material is aluminum, 


with a thermal conductivity of ky = 190 w/m—° K. 


Popsit 
SEST imTEAPRE TATION 


| 
O (4m G/u) | | 
1$ 20 30 40 SO 6070 





Fin pitch = 2.0 per in = 78.74 per m 

Plate spacing, Db = 0.750 in = 19.05 x 10°3m 

Fin length = 12.0 in = 304.8 x 10°3m 

Flow passage hydraulic diameter, 4r, = 0.0474 ft = 14.453 x 10°3m 

Fin metal thickness = 0.032 in, aluminum = 0.813 x 10°3m 

Total heat transfer area/volume between plates, 6 = 76.1 ft?/ft? = 249.672 m?/m?3 
Fin area/total area = 0.606 


Plain Plate-Fin Surface 2.0 
Figure 23 
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Chapter 3 


CONVENTIONAL 
DESIGN METHOD 


3.1 Introduction 


This chapter outlines the conventional design procedure that was used to 
size the single-pass, counterflow gas-to-gas heat exchanger. This procedure 
was performed for various combinations of heat exchanger effectiveness and 
specified pressure drops. A complete numerical example is provided in Ap- 
pendix B. Final results for this conventional method are shown in Table 
Sal. 


3.2 Conventional Methodology 


For a single-pass counterflow design, when the core dimensions on one side 
are fixed, the dimensions on the other side are also known. This means the 
design is driven by the side that has the more stringent AP/P requirement. 
The method of determining the controlling side is contained in Appendix 
B. In every case considered, the hot side of the heat exchanger was the 
controlling side and so its requirements drove the design. 

The following is a step-by-step heat exchanger sizing design procedure 


which closely follows that outlined in reference [2]. 


1. For the given heat exchanger effectiveness, determine the fluid outlet 
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temperatures. Calculate the fluid mean temperature on each side and 


evaluate fluid physical properties p;, po, (1/p)m, and uw. 


2. Determine N;,, for the exchanger and then N;,, for each side. The in- 
fluence of longitudinal heat conduction is ignored in this first iteration 


of design. 


3. Estimate hot and cold side pressure drops, select an appropriate Np,z, 


and then a value of 7/f from Figure 2.3. 


4. Calculate mass velocity, G, from information in steps 1-3 and the cor- 


responding value of AP/P. 
5. Calculate Np, and determine values of j and f from Figure 2.3. 


6. Compute heat transfer coefficient, h; temperature effectiveness of the 


fins, n¢; and the total surface temperature effectiveness, no. 


7. Calculate heat transfer area, A; minimum free flow area, A.; heat ex- 


changer frontal area, Ay+,; flow length, L; and heat exchanger volume, 
Vol. 


8. Compute AP/P from known conditions and calculated parameters. 


3.3 Analysis 


For each case show in Table 3.1, cycle efficiency with losses, 7, was calculated 
and the cycle efficiency degrade, no. — 7, is listed for each case. A correlation 
between heat exchanger effectiveness, €, total allowed pressure drop, k, cycle 
efficiency degrade, and heat exchanger volume was sought. 

The relationship between cycle efficiency degrade and effectiveness at 
various levels of k is shown in Figure 3.1. The increased level of degrade 
at lower values of effectiveness is as expected, and for a particular value of 
effectiveness, the degrade increases with increasing pressure drops. 

The effects of varying effectiveness and pressure drops on heat exchanger 


volume is depicted in Figure 3.2. The higher the effectiveness, the longer 
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and therefore, larger the heat exchanger. Similarly, the greater pressure 
drop for a particular effectiveness also increases the volume. 

Figure 3.3 was derived by extracting values of effectiveness at nominal 
values of 799 — 7 and plotting them on the appropriate k curve of Figure 3.2. 
This resulted in curves of constant 7~ — 7 that have a minimum volume 
at a particular value of k and e. It is emphasized that these “minimum” 
volumes were obtained at constant values of cycle degrade with varying 
pressure drops and effectiveness. These “minimum” volumes are plotted 
against effectiveness in Figure 3.4 with the respective k values annotated on 


the graph. 
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Table 3.1: Conventional Method Results 
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Efficiency Degrade vs Effectiveness 
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Chapter 4 


MINIMUM ENTROPY 
DESIGN METHOD 


4.1 Introduction 


The use of irreversibility concepts in the design of the gas-to-gas counterflow 
heat exchanger will be presented in this chapter. For the same mass veloci- 
ties that were used in Chapter 3, the heat exchanger geometric parameters 
can be determined to arrive at minimum levels of entropy generation. The 
outputs of this design procedure are heat exchanger effectiveness and the 
pressure loss suffered on each side. For a more detailed explanation of this 


procedure, the reader is directed to references [4] and [6]. 


4.2 Background 


Since the method deals with the entropy generation within the recuperator, 
it is appropriate to begin with an entropy flux analysis on the control volume 
of the recuperator in Figure 4.1. The heat transfer from the outer walls is 
assumed to be negligible. The entropy generation rate can then be written 
as 


S= m1(Sout — Sin) + M2{Sout — Sin)e (4.1) 
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Stream 1 


Length 


Ta 
cn 


Recuperator Control Volume 


Stream 2 


Figure 4,1 


Substituting expressions for the ideal gas entropy changes based on pressure 
and temperature relationships, this becomes 


Ss = Cr anltalitoce tan) F (R/Cp)1!n( Pi/Prout)| 


+Cmaz|!n(T2oue/T2) = (R/Cp)eln( Po / Prout) | (4.2) 
The number of entropy production units, N,, is defined as 
Np= S/O (4.3) 


Writing expressions for the first law of thermodynamics of the control vol- 
ume, 


Crone i Trout) te Crane ca Teaut) = 0 (4.4) 


and heat exchanger effectiveness, 


€ = (isu 7 T,)/(T2 a Ti) (4.5) 
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permits the elimination of Tyou: and T2ou¢ in equation (4.2) which allows the 


rate of entropy production formula to be nondimensionalized as 


C'mi ly Cmin i 
N, = Sin [1 +e (2 —1)] +i [1 - (1-2)| 


maz 1] maz T2 


Cue (4) ( = (Z) a ( — (4.6) 
—- 71 —_ — —= tre: c 
Cras Nae Pome Cs) Poe 


A special form of this last expression can be obtained in the case of nearly 
ideal heat exchangers. When the stream-to-stream AT’s and frictional A P’s 


are held small, two inequalities can be applied. 








l-e <<1 (4.7) 


(AP/P)12 <<1 (4.8) 


Applying these and expressing € in terms of Nt, yields 


Omi n T2 Cmin ( Tei ) 
Cmaz Gi @ nae T» 


I 








+ (See) (1 - Gain) (1- R) exp |-Nuu (1- ga) | 
Cmaz C'maz 1 — Gain. (1-H) 1 - Sminexp |—N,, (1 — Saia)] 


Cr. R IE R AP 
- Gain (B) (AP) , (8) (42) on 
The first two terms can be characterized as a contribution from capacity 


rate imbalance, the third represents contribution due to finite N;,,, and the 


last two show fluid friction effects. 





Requiring the streams to be balanced and applying the calculus of limits 


as Cmin — Cmaz, equation (4.9) can be reduced to 


Cmin l T2 


2 
~~ in— + in [1 - Sm (1-2) 4 a2 i oe 
Cmaz T Cmaz T> dA T>2 N tu 
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sue (2) (3) +(2) (HF), wo 


When the overall Ni, is expressed in terms of N;,, for each side of the heat 


transfer surface, 
1 1 ONE | 


Ntu - Neu Cmaz Ntu2 


the number of entropy generation units can be divided into three contribu- 








(4.11) 


tions 
N, = fe imbalance + N51 of N32 (4.12) 


Applying the limit Cmin — Cmaz, the imbalance terms of equation (4.10) 
can be written as 


Cmaz T: Car T: 
iy Newnbalance == ( a 1) & —-i- G n= | (4.13) 








which will clearly vanish when Cmin = Cmaz. Since this is the case for the 
design here, the number of entropy production units for each side can be 


determined from 


2 
Es abs 1 R iS 
Nee Z ae at 4.14 
es ( ence ) Nea’ (=) Bee ate) 
IT. [T. ) 1 R AP 
N, = 2 1 eee (=) a5 
z ( 7B Ne (=) P /, ely 


4.3 Design Considerations 


The number of entropy generation units that have been deemed acceptable 
for a particular design is a function of the hydraulic radius, minimum free 
flow area and flow length. It is necessary, then, to express N, in terms of 
these and other known flow parameters. Using the definitions of Ni, and 


friction factor for each side 


Ntui,2 (L/rn) Not (4.16) 


(AP/P)1,2 f(L/rn)G*/(2pP) (4.17) 
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in equations (4.14) and (4.15) leads to 


ar 
Ne ee EO TL ; 4, 
1,2 ae f{(L/rn)g (4.18) 
where 
ay= I) 
a2 = Oe Cmar (4.19) 
by a Cn Caras 
by = 1 (4.20) 
By = (R/C,): 
3 (R/C5)2 (4.21) 
T: a ° 
= me Se 
—— ( 7 \/ 7, | (4.22) 
I= G/(2pP)2 (4°23) 


As can be seen in equation (4.18), the ratio L/r, performs a trade-off func- 
tion, i.e., for a fixed g and Npr,, there will be an optimum L/r;, which results 
in a minimum N,;2. When the function Nzi,2 is minimized, the (L/rpj) opt 
is given by 


@ Z T/Nst z 
(L/rrJopt = G (sre ri Ta) ee) 


4.4 Minimum Entropy Design Methodology 


The method used is based on optimizing (L/rp,). Since rp is set by the selec- 
tion of the plate-fin surface, this yields a length that will produce minimum 


entropy for the given hydraulic radius. 


27 





The following outlines the methodology employed. A complete numerical 


example is shown in Appendix C. Final results for all cases are tabulated in 


Table 4.1. 


1. For the cases listed in Table 3.1 that resulted in a heat exchanger 


volume, use calculated G and Np, to determine (L/rp) opt. 


2. Compute resultant heat exchanger effectiveness and hot and cold side 


pressure drops. 


3. Calculate N;,, Ns, and 


Stortal’ 


q 


4. Calculate the heat exchanger volume for each case. 


5. Determine the cycle efficiency degrade, No — 7. 


4.5 Analysis 


Figure 4.2 shows the nearly linear relationship of cycle efficiency degrade and 
heat exchanger effectiveness. There is not the clear definition of different 
values of pressure drop as was observed for the conventional method in 
Figure 3.1. In fact, there is even a wider range of k values here, from 
k = 0.056 for case E-6 to k = 0.269 for case E-24, and yet the results are 
almost linear across that range. 

The volume plotted in Figure 4.3 is that volume obtained by minimizing 
entropy generation in the heat exchanger and should not be interpreted as 
a minimum volume, but rather as an “optimum” volume for minimizing 
entropy generation. 

To achieve a better basis for comparison of the two methods, the entropy 


generated by the various designs must be discussed. 
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| Case # Lopi(m) Vol(m*) k € Noo — 7 Nea 


E-6 

E-7 

E-8 

E-9 

E-11 
E-12 
E-13 
E-14 
E-16 
E-17 
E-18 
E-19 
E-21 
K-22 
K-23 
E-24 


32.03 
25.86 
22.04 
20.12 
24.35 
19.08 
16.11 
15.08 
19.80 
15.40 
13.19 
12.11 
16.98 
4 feiel 
11.52 
10.47 


0.0563 
0.0727 
0.0876 
0.1014 
0.07658 
0.1069 
0.1320 
0.1559 
0.996 
0.1422 
0.1778 
0.2058 
0.1252 
0.1799 
0.2251 
0.2685 





Table 4.1: Entropy Generation Method Results 
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Chapter 5 


COMPARISONS AND 
RESULTS 


5.1 Introduction 


To this point, the minimum entropy method results have been compared to 
the conventional method results in terms of conventional method standards, 
1.e., pressure drops, efficiencies, and effectiveness. It is necessary to compare 
these methods in terms of the minimum entropy method, 1.e., the entropy 


generated in each design. 


5.2 Conventional Method Entropy Generation 


The number of entropy generation units for each design can be computed 
from previously determined information and equations (4.14) and (4.15). 
An example of these calculations is contained in Appendix D. Results for 
the conventional cases are listed in Table 5.1. The large deviation in N, and 
Noo — 7 for cases C-21 through C-24 are a result of violating the assumption 
that 1 — ¢€ << 1. These conventional cases will be discarded from further 
consideration. However, the minimum entropy cases at corresponding mass 
velocities are still valid. 

When the entropy generation units are plotted against the degrade in 


cycle efficiency, Figure 5.1, the dependence on pressure drop is still present. 
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The grouping of data points is caused by the discrete changes in heat ex- 
changer effectiveness, decreasing from left to right. 

Figure 5.2 was created by using the same values of 7,, —7 that were used 
for Figure 3.4. The value of N, at each of these cycle efficiency degrades 
was read off at a particular value of k. Plotting these “minimum” volumes 
against entropy generation rate exhibits the general trend that the larger 
the volume of the heat exchanger, the lower the total entropy generated. 


This is in consonance with the results predicted by reference [4]. 


5.3 Minimum Entropy Generation 


The cycle efficiency degrade vs. entropy generated for the minimum entropy 
method 1s similar in overall shape to that of the conventional method, but, 
as Figure 5.3 shows, the relationship is nearly linear and does not show the 
discrete pressure drop differences. 

Figure 5.4 shows that the “optimum” volumes determined by(L/rh) ope 
follow essentially the same path as that of the conventional method. For a 
direct comparison, Figures 5.2 and 5.4 are plotted together on Figure 5.5. 

To ensure that the (L/rp)ope expression, (4.24), was giving minimum 
entropy results, a series of calculations was performed at (L/r;,) values on 
both sides of the optimum for case E-9. Figure 5.6 shows that the minimum 
entropy solution was, in fact, being determined. It is also noted that the 
magnitudes of N, and 7. — 7 are essentially the same over this range of 


heat exchanger effectiveness. 


5.4 Comparisons 


The calculations performed on case E-9 to prove minimization also yielded 
heat exchanger volumes at a constant mass velocity G. These additional 
data points are plotted on Figure 5.7. The dash-dot line connecting them 
is a locus of recuperator volumes along a constant G path. This shows why 
the volumes at minimum N, or 7. — 7 along a constant G path are larger 


that the miminum volume of a given N, or Noo — 77. 
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Figure 5.8 shows a comparison between conventional “minimum” vol- 
umes and minimum entropy “optimum” volumes. The minimum entropy 
method yields volumes that are larger than the conventional method for the 


same cycle efficiency degrade. 
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Case # Non 





Table 5.1: Conventional Method Entropy Generation 
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Chapter 6 


CONCLUSIONS 


6.1 Discussion 


The nature of the results and conclusions drawn for the example cycle would 
be the same for other cycles and recuperator surfaces. For high effective- 
hess recuperators, the counterflow arrangement is dictated. Here with the 
same finned surface on both sides, the (AP/P)./(AP/P)a = (P,/P.)* or 
approximately 1/25. If a closer plate spacing were used on the clean air cold 
side more AP/P would appear on the cold side for any total UP/P. This 


would result in even smaller minimum volumes at any given N, or 7 — 7. 


6.2 Conclusions 


In selecting a gas turbine recuperator, essentially the same results are ob- 
tained either by the conventional method of minimizing degradation of cycle 
efficiency, Nog — 7, or by minimizing entropy generation, Ns. 

The selection of an “optimum” recuperator requires that care be exer- 
cised in the selection of the path along which the optimization is calculated. 
It has been shown that the volume of the recuperator for minimum JN, or 
No — 7 along a constant mass velocity path is much greater than the true 


minimum volume at any No — 7 or Ng. 
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Appendix A 


Sample Calculation for 
Effect of Recuperator on 
Cycle Efficiency 





1 
eee 

n Wnet Wt + We 
Co = : — : 
Q Q 


i = m(hg — hs) = me,(Ts — Ts) 
We = m(hy = he) = mcep(Ti a T2) 


Q = m(h4 — hg) = mne,(Ts — T3) 
Turbine: 


T4 (Bayo) 


Te oe 


es 
Ps (npe)( 2+) i 
_ o | 


T4 
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(.9)(444)] 7 
ip [a = 820.1°k 
Compressor: 
T> _ (=) EG 
71  £4=\PR 
-1 
P, Tipe 
i — 
2 ai (3 
T, = (300)(6) 0-7 = 514.2°k 
Recuperator: 
ee 
— Ts — Ty. 
i (e) (Ts _ T2) + T5 
Tz; = (.8)(820.1 — 514.2) + 514.2 = 695.2°k 


Cycle Efficiency: 


ein Foe 1200 1820-1) -- (300 — 514 2) 
0 
is 1300 — 695.2 
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Appendix B 


Conventional Design 
Method Numerical Example 


Before calculations can be performed, two points need to be clarified: 1) the 
method of determining how loss of turbine work due to pressure drops in 
the system is distributed and; 2) the method for determining which side of 
the heat exchanger has the more stringent AP/P requirement. 


Turbine work is given by the expression 


P oy 
fameps [1 ~ (FZ) | 


we 


P 
eee 1 (22) Weer 
mnicpTh P,—-AP,/ Ph} }1— 95h 








Expanding the last term, 

i cae NP AP. 7 
poet ~ 1- (=) fi + S44 Ste] 
MmNtCplp Ten Py i. 

aN (Ny ag AP. 
MnNtCpTp Ph Ph P, 


so that loss of turbine work can be expressed as 


Aw oe (F:) (44+ 4) 
mmCpTh ap eR iz 
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and it is clear from Figure B.1 that oor determines the loss of turbine 


work irregardless of how it is distributed. Introducing the definitions r = oh 


TEMPERATURE 





specific entropy 


Temperature —- Entropy Diagram 
Figure Bl 


and 9 = (#) , 1t can be shown that for a counterflow exchanger where 
¢/ av 
m= Ap, and L, = Lp, 








ia eed 
a) op 
Using a slightly different form of the last turbine work expression, 
Aw: _ 1 oe ri 1 I) Ef 
mmcpT, rt LAP,/P, P, 
eee (441) 2 
mncpTh srt \ P. 
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: Ae AE 
Introducing apes PL 


drops are, 


Awt 


We 





Aw: 


We 





(6/r? +1) AP, 
(nt) PF. 





1 
es) 


(0/r? +1) AP, 
r?—1 £P, 





= k, expressions for hot and cold side pressure 


AP, _ k(6/r?) 
Ph 1+ 6/r? 
WE k 
Ps 14+ 6/r? 


The numbered steps that follow correspond to those outlined in Section 


3.2. The input design parameters used for the example calculations were 
e— 0-8 and & = 0.04. 


B.1 Temperatures and Fluid Properties 


so that 


i = 


R= 


Ps 


T, = 300°k 
i a 
Tins Peer) 


E — 
Ti | — 
P, 


1.41 
300(5) G5)U-47 = 486.8°k 


P. 
—(1—k) =5(1— .04) = 4.8 
iat 
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Ts © 


Js 


Ts = 





1.4—1 —1 
(4.8)¢°)(AS 
(i ————_____--—_ = 868.5°k 
: | 1300 
os Ts — Te 
Si 7; 
so that 
T3 = e(Ts = T2) ate T» 
Tz; = (.8)(868.5 — 486.8) = 792.2°k 
Similarly 


T6 = Ts as (Ts 7 T2) 


Ts = 868.5 — (.8)(868.5 — 486.8) = 563.1°k 


Since this is a counterflow design with C* > 0.5, the best choice for 


average temperatures is the arithmetic average on each side 


Ts+Ts 868.5 + 563.1 


it ad — (liga: 
hm 9 9 
486.8 + 792.2 
Icom = ies tts = gece (ee a = 009-5 0k 
2 2 
Inlet, outlet and average densities are next computed. Pv = RT or 


ae, © butv — 1/p,soig— P/ iia 


1.126 x 10°N/m? 


PS ~  (287j/kg —° k)(868.5°k) 


= 0.4517kg/m? 


90 





1.126 x 10° 
—- 7 =06967k¢/m? 
oe (287) (563.1) g/m 


1-7-1 1 
(=) = = (— + ~) = 1.824m?/kg 
Dyan 2\ps pe 


5.63 x 10° 
i ee A O8Oka/m* 
a (287) (486.8) g/m 


5.63 x 10° 
Se Re 
i. (287) (792.2) g/m 


; . 
@ = 0.3260m°/kg 
P/ mye 


’ 


ee LV 


For a couterflow heat exchanger, 
eerie c- | 
~ 1—c*ezp[—Ni,(1 — c*)] 
but for the balanced flow case (c* = 1) under design, 
Neu 
i 
1 ate New 





or 


Since gas is the working fluid on both sides, it was estimated that both sides 


would have approximately the same surface resistance so that 


Neu h = Neue =~ 2Niu = 8 


B.3 Estimated Pressure Drops 
Recall 


Tyee, 168015 


ol 





and 





ce = a 27.983 
fal 2) 
AP; ES AC ee = (0.04)(27.983) — 0.03862 
ee 1+ 0 /r? 1 + 27.983 
AP. k 0.04 — 0.00138 


P,  1+06/r? 1+ 27.983 


From Figure 2.3, an estimate of j/f = 0.0035/0.0077 was obtained. 


B.4 Mass Velocities 


The mass velocity on each side is estimated from the relationship 


( 29-Pre )) i) 
(Up) aes Neu stde if 


As a first approximation, 7, was assumed to be 0.90. 


G, = (ee ee) — a, 


Gaz 

















(1.824m3/kg)(0.7 8 0.0077 
Gy, = 17.586kg/m? — s 
1 
— (es x 109) = (208) 2 
(.3260)(0.7)3 8 0.0077 
G. = 17.583kg/m? — s 
Fy _ GrDp _ (17.583kg/m? — s)(0.0474f t)(.3048m/ ft) 
oe Lh 7 3.396 x 10-°kg/m —§ 
NRe,h — 7481 


O2 





(17.583)(.0474)(.3048) 
3.159 x 10-5 


Entering Figure 2.3 with these values yields 


NRe,c = OU4Z 


fr = 0.0078 
jn = 0.00355 
fe = 0.0077 
je = 0.0035 


B.5 Heat Transfer Coefficients and Fin Effective- 
nesses 


With the Colburn factor known, the heat transfer coefficient can be deter- 
mined from 


= 2 
h = jGepP, ° 


This result is used to compute the fin characteristic length, 


L 
(2h)? 
 \ 6 


The temperature effectiveness of the fin is next calculated from 


tanhml 
n= 





ml 
A good approximation for / is given by half the plate spacing minus the fin 
thickness. 


The total surface temperature effectiveness can then be calculated 
Mo =1—(1—nz)Az/A 
where A¥f/A is a characteristic of the surface selected. 


(0.00355) (17.586kg/m? — s)(1.0k7/kg —° k) 


hp 5 
(0.7)3 





Ah 


M< 


"fc 


"o,c 


79.19w/m? ~°k 


(2)(79.19w/m? —° k) 2 
(190w/m —° k)(0.032in)(0.0254m/in) 


32.024m7! 





cf 
€ an _ 032in | 0.0254m/in 


8.71 x 107*m 


tanh(32.024)(0.71 x 1073) 
(32.024)(0.71 x 10-3) 


0.9749 


1 — (1 — 0.9749)(.606) 


0.9848 


(0.00355) (17.583)(1.0) _ 7. 9. Im? k 


(0.7) 3 
L 
ee ee siren! 
(190)(0.032) (0.0254) 
~3 
tanh(31.79)(8.71 x 107°) ce 


(31.79)(8.71 x 10-3) 


1 ~ (1 — .9752)(.606) = 0.9850 


The overall heat transfer coefficient can now be determined neglecting foul- 


o4 





ing resistances and wall thermal resistances, 








1 . en! 
U (noh)n  (noh)e 

1 . 1 1 

= = ee 0.02583 
u (0.9848)(79.19) ' (.9850)(78.06) 


U =38.714w/m* -°k 


B.6 Dimensions 


The heat transfer surface area is computed from 





Nin? 
A= 
ue 
Ok7/kq —°k 
A = (AN(8600kg/hr)(1.0k9/kg -* ) _ 59 orn? 


38.714w/m? —°k 
Minimum free flow area is given by 


A, = (W/G) = —————__ = 0.05686 
US) 17.586kg/m? — s i 
In order to determine frontal area, a must be computed from the relationship 
— dBD;/4 
— d+d4+2t 


Here a parting plate thickness of 0.4mm was assumed, so 


(0.75in)(76.1f t? / ft3)(0.0474/4 ft) 


7 0.75in + 0.75in + 2(.4mm)(.1/2.54 in/mm) 
Go =F 0.44162 
Afr = A./a 
0.05686m? 
A = ———-—O0.1 ? 
fr 0.44162 aon 
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The flow length is given by 


_ [DA 

a Ga 

z (0.0474 ft)(0.3048m/ ft)(103.321m7?) 
7 (4)(0.05686m?2) 

L = 6.56m 


The total volume of the heat exchanger is calculated by 


Vol = (Ay,)(Z) 


Vol 


(0.12876m?)(6.56m) = 0.845m* 


B.7 Calculated Pressure Drops 


For this level of design, entrance and exit effects are not considered, so only 


core friction effects will be used. 





I 
> = Hats aa) rer 2 
AP,, (0.0078) (6.56)(1.824) (17.586)? 
SS = 0.00125 
P,, (2)(0.0474/4)(.3048)(5.63 x 105) 
Similarly 
ie 0.00125 
so that 


Psi 
4 —— = 0.03674 
P 
The cycle efficiency with these input parameters is 


(T, — Ts) + (T, — To) 


! 14 —-Ts3 


(1300 — 868.5) + (300 — 486.8) 
(1300 — 792.2) 


06 





n = 0.4819 
and the efficiency degrade is given by no — 7. 


Noo — 7 == 0.5760 — 0.4819 = 0.0941 
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Appendix C 


Minimum Entropy Design 
Method Numerical Example 


This example uses the data from case C-12 of Table 3.1, 
G = 12.456kg/m? — s, 7 = 0.00355, and f = 0.0078. 


C.1 Optimum Length 


In order to use equation (4.24), intermediate values must first be determined 


a 0.00355 
Nse= oso = os =: 0.00450 
= We (O77 
The temperature span parameter, 7, cannot be calculated until Ts is known, 
so an estimate based on previous iterations is used. If the calculated value 


of r is close to the value used, the computations are valid. 


1 
2 r/ Not 2 
a 
Ul ra)ont = GBF] peP. + 1] onPh) 
2 0.365 /0.00450 : 
(L/rr)opt = 12.456 ! 0.325 1.818 
(0.287) (0.0078) pe T Ti7ex — 


(L/ra)opt = 5283 
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Since r, is set by the geometry selected, L can be calculated 


Lopt = 5283(.0474f t)(.3048m/ ft) 
opt = 19.08m™ 


C.2 Heat Exchanger Effectiveness and Pressure 
Drops 


From equation (4.16), 


Nish = (L/ta) Nee 


Niun = (5283)(0.00450) = 23.79 
Applying the same assumption as was used in Chapter 3, 


Ntuh ae Niue = 2Neu 


New 23.79/2 = 11.89 


Then, since this 1s still a balanced counterflow design, 


For hot and cold side pressure drops, equation (4.17) is used. 





AP 

(=), = f(L/rn)9°/(2eP) = f(L/rn)9” 

AP ; 

(=-)h = (0.0078)(5283)(0.05005)* = 0.010322 
AP , 

(=5-)e = (0.0077)(5283)(0.00946)’ = 0.00369 
AP 


og 





C.3 Entropy Generation 


Equations (4.14) and (4.15) will be used to determine entropy generation on 


both sides. 
r? R AP 
rw = (8), (8) 
: Ntuh (2) 1s h 


0.365 
a: 103 
Nsp 55-79 + (0.287)(0.10322) 





Ns, = 0.04497 


07365 
N = — 0.287)(0.00369 


Ns, = 0.0164 
Ns = 0.04497 + 0.0164 = 0.0614 


C.4 Heat Exchanger Volume 


Since Af, is calculated from minimum free flow area and other surface geo- 
metrical constraints,Ar, from the corresponding case in Chapter 3 was used 


to calculate the heat exchanger volume. 


Vol = L(Afs,) = 19.08m(0.12876m?7) = 2.46m° 


C.5 Cycle Efficiency Degrade 


Temperatures that reflect the actual effectiveness and pressure drop must 
first be determined. 
ii P, 
— —(1—k) = 5(1 — 0.1069) = 4.465 
B= FUR =5( ) 
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(5 le 
i Ps 

5 i 

(4 465) (o\(ieseaye ) 
1.4 

TT = ee = 884.8°k 

: | 1300 | 
T3 = «(Ts = T2) ai T» 
is) = 0.922(884.8 _ 486.8) + 486.8 


a= 1Soekek 


(T, — Ts) + (T1 — To) 


ia Ty oe 
(1300 — 884.8) + (300 — 486.8) 
Ee eee ee) =. 0) 5119 
7 1300 — 853.8 
ee oe GOO or 006412 


A check on the validity of the r value assumed shows 


2 
a Te aR 
¥ T) Ts 
z 
— [884.8 [486.8 
7 486.8 884.8 


0.367 


~~ 
i 


Therefore, assumed value of 0.365 was valid 
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Appendix D 


Sample Calculation of 
Entropy Generation for 
Conventional Designs 


Example uses case C-12. From previous example, 


rT = 0.34461 
AE, 
— =_ 0.0359 
ef 
Ae 
—— = Q.00136 
ie 


Niuh aa Niue = 8 
Using equations (4.14) and (4.15), 





N _ T a R (>) 
3,h Ntuh C A P h 
0.3446 
Nzn = eeuaee a (0.287)(0.0359) = 0.05338 
Similarly, 
0.34461 
Nsc = ia + (0.287)(0.00136) = 0.04347 
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N, = NenpatNs,- = 0.09685 
Adding the first term of each side’s expression gives the AT contribution, 


and the sum of the second terms is the AP contribution to the entropy 
generation. 


Neat = 0.08615 


N;.aP 0.01070 
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